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Abstract. We present a new analysis of the very first mid-infrared N-band long-baseline interferometric observations of an 
extragalactic source: the nucleus of the Seyfert 2 galaxy NGC 1068, obtained with MIDI (Mid-InfrareD Interferometer), the 
mid-infrared beamcombiner at the European Southern Observatory (ESO) Very Large Telescope Interferometer (VLTI). The 
resolution of T/B = 10 mas allows to study the compact central core of the galaxy between 8 and 13 pm. Both visibility 
measurements and MIDI spectrum are well reproduced with a simple radiative transfer model with two concentric spherical 
components. The derived angular sizes and temperatures are about 35 and 83 mas, and 361 K and 226 K for these two com¬ 
ponents respectively. Some other evidences strongly support such low temperatures. This modeling also provides the variation 
of optical depth as a function of wavelength for the extended component across the N-band pointing towards the presence of 
amorphous silicate grains. This shows that MIDI actually carried out the first direct observations of the distribution of dust 
around the central engine. Together with other observational pieces, we are able to move a step forward in the reconstruction of 
the picture drawn for AGNs. 
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1. Introduction 

NGC 1068 is one of the brightest and nearest Seyfert 2 galax¬ 
ies. Located at a distance of 14.4 Mpc, it is unique for the study 
of the active galactic nucleus (AGN) it harbours. The nucleus 
of NGC 1068 has been well studied in the entire spectrum, 
from X-rays and UV to radio wavelengths, including the op¬ 
tical, near- and mid-infrared wavelength ranges. The location 
of all the components emitting at different wavelengths has 
been revisited in Galliano et al. <f2003l> . NGC 1068 was clas¬ 
sified as a Seyfert 2 with respect to the narrow emission lines 
it emits. However, Antonucci & Miller t 19851 have discovered 
broad, polarized emission lines, which suggest the presence of 
a Seyfert 1 nucleus, hidden by a geometrically and optically 
thick dusty torus surrounding it. Thus, they laid the foundation 
stone of the unified scheme of AGN. Several theoretical stud¬ 
ies of AGN tori were reported by Pier & Krolik ( 119921119931 . 
Granato & Danese IMS. Efstathiou et al. (TT9951 . All these 
models consider radiative transfer in an anisotropic, but ho¬ 
mogenous dusty torus. Nenkova et al. 12002 1 considered radia¬ 
tive transfer between several clouds of dust along radial rays 
through the torus, and overcame some difficulties of the former 
homogeneous models for the description of the spectral energy 
distributions. More recently, Schartmann et al. d20r)5l> describe 
the three dimensional treatment of transfer radiative through 


dusty tori at hydrostatic equilibrium and succeed in reproduc¬ 
ing both the mean large aperture spectra from UV to far-IR of 
several Seyfert type 1 galaxies, and the recent MIDI observa¬ 
tions of two specified Seyfert type 2: Circinus and NGC 1068, 
except for the 9.7 pm silicate feature. Geometrical and dy¬ 
namical properties of dusty tori are theoretically investigated, 
among others, by Krolik & Begelman ( 119881 1. Zier & Biermann 
(Emu, and Beckert & Duschl ( 2003 ). 

An alternative explanation for the unified scheme of AGNs 
is given by Elvis ( 2000 .). who actually proposes a simple 
and empirical scenario of funnel-shaped thin shell, warm and 
highly ionised outflow that could account for all the struc¬ 
tures in the inner regions of AGNs with various luminosities. 
However, this model is not in contradiction with the obscuring 
torus paradigm since the accelerating bi-conical wind could be 
itself considered as a form of dusty torus. 

In the near-IR, the first J, H and K-band Adaptive Optics 
(AO) images of NGC 1068, at 0.12” resolution (1” corresponds 
to 72 pc), were reported by Rouan et al. ( TPMl l. They gave 
an upper limit on the size of the unresolved K-band core of 
~ 9 pc, and highlighted an elongated S -shaped structure ra¬ 
dially extending up to 20 pc, at a position angle (PA) ~ 102°. 
Spectroscopic high angular resolution X-band observations ob¬ 
tained with PUEO-GriF at CFHT are presented by Gratadour 
et al. ( 2003 ). Emission of hot dust seems to be confined in an 
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almost resolved region of 120 mas. Based on K, L, and M-band 
diffraction-limited images of NGC 1068 obtained with NACO 
at the VLT, Rouan et al. ( 2004 1 reported several emissions 
around the nucleus being structured in several regions. The K 
and L-band emissions are actually resolved with a FWHM of 
67 and 122 mas respectively, and are elongated in the NS di¬ 
rection. 

First 76 mas resolution K -band images of NGC 1068 
obtained from bispectrum speckle interferometry were inter¬ 
preted in terms of a resolved gaussian component of FWHM 
~ 30 mas (Wittkowski et al. ll998t >. Speckle imaging performed 
by Weinberger et al. (il999t at the Keck Observatory led to 
agreeing results. 

Then, the first K-band long baseline interferometric mea¬ 
surement of NGC 1068 was obtained with VINCI at the VLTI 
(Wittkowski et al. 120041 . Coupled with K-band speckle inter¬ 
ferometry measurements (Weigelt et al. 120041 . they account for 
the observations with a two-component model where the small 
component has an angular size of less than 5 mas and the other 
one has a size of 40 mas. Weigelt et al. (120041 present near-IR 
bispectrum speckle interferometry made with the S AO 6m tele¬ 
scope. The K'-band FWHM diameter of the observed compact 
core is ~ 18x39 mas, and the PA of the north-western elonga¬ 
tion is -16°. In the //-band, the FWHM diameter of this same 
component is ~ 18x45 mas, and the PA is ~ -18°. 

First mid-IR images with 0.1” spatial resolution are re¬ 
ported by Bock et al. ( 2000 ,). using the 10 m Keck II telescope. 
They interpret the northern elongation of the central peak as 
reprocessed radiation from the AGN. Tomono et al. ( 2001 1 
also obtained mid-IR images with 0.1” spatial resolution with 
MIRTOS on the 8.2 m Subaru Telescope and confirm that the 
central region of NGC 1068 is elongated in the NS direction, 
with a FWHM of 0.29” x 0.18”. They use a modified grey-body 
radiation, with absorbing silicate features in order to model the 
spectral energy distribution (SED) of the central peak. 

Jaffe et al. ( 20044 report the first mid-IR /V-barul long- 
baseline interferometric measurements using MIDI, the VLTI 
mid-IR beamcombiner. They consider a model composed of 
two ellipsoidal gaussian disks. They derive sizes of about 
10 + 2 mas in the NS direction, parallel to the inner ra¬ 
dio jet, and <12 mas in the direction transverse to the jet 
for an inner hot component (T > 800K) and sizes of about 
30 + 5 mas x 49 ± 4 mas for an outer warm component 
(T ~ 350 K). 

Infrared interferometry clearly has the potential to study the 
circum-nuclear environment of AGNs to bring evidences for or 
against the unified sheme to help better understand their nature. 
In this paper, we report an alternative analysis of the first N- 
band MIDI observations of the nucleus of NGC 1068. First, 
in Sect. |2] we give a short review of the observations and data 
processing. Sect.0details the models considered in this paper. 
Then, the results are presented in Sect.0] and further discussed 
in Sect. 0 Sect. Ogives the general conclusions of the study. 

2. Observations and data reduction 

The observations were performed in 2003 (June 15th, 
16th; November 9th, 10th), during four nights of Science 


Demonstration Time with MIDI, the mid-IR /V-band (8- 
13 /jm) instrument. Observations were carried out with two 
8.2 m diameter Unit Telescope (UT) at the European Southern 
Observatory (ESO) VLTI on Cerro Paranal. A 0.6” x 2” slit was 
used with a low resolution prism, set in the direction perpendic¬ 
ular to the slit, providing a spectral resolution of A/ AT ~ 30. 

MIDI has two interferometric outputs which are the sum of 
the flux of NGC 1068 and the fringe modulation, plus the back¬ 
ground. The NGC 1068 flux and the modulation were measured 
separately as the photometric outputs were not operational at 
the time of the observation. The two interferograms are com¬ 
plementary with opposite phases. Using this property, the two 
outputs are subtracted to disentangle the modulation from the 
continuous signals. The amount of continuous flux from NGC 
1068 was accurately obtained thanks to the “chopping” tech¬ 
nique: both UTs alternatively point to the sky and to the source. 
Subtraction of the sky sequences from the source sequences 
yields the source flux. The chopping frequency was 2 Hz dur¬ 
ing the two nights. Since no spatial filter was used at that time, 
the flux from NGC 1068 can be considered constant. The mod¬ 
ulated part of interferograms is then normalized by the contin¬ 
uous flux. The contrast of individual fringe scans is obtained 
with a classical method which has been used with precursors 
of VLTI (Coude Du Foresto et al. 1 997 1 . Normalized scan se¬ 
quences are Fourier transformed. The power spectrum is band 
pass filtered in each MIDI spectral channel. The squared fringe 
contrast is the integration of the filtered power spectrum from 
which the photon and read-out noise bias have been subtracted. 
Statistical errors on squared fringe contrast estimators are set 
from the dispersion of measurements. 

2.1. Visibilities calibration 

The entire data reduction has been done with the software de¬ 
veloped at LESIA at Paris Observatory for MIDI 1 . It also re¬ 
lies on methods developped for precursors of the VLTI (Perrin 
l2003il . This software has been used to reduce data on stars 
for which visibility functions can be more easily anticipated 
(Kervella et al. 12005 ? accepted). 

As the pupils of the UTs are larger than the spatial coher¬ 
ence diameter of the wavefront, and as there was no spatial fil¬ 
tering at the time of the observations, beams are strongly sen¬ 
sitive to turbulence and alignment defects, causing coherence 
losses. Therefore, visibilities need to be carefully calibrated. 
To do so, a calibration star is observed before and after the ob¬ 
servation of the source. In the case of NGC 1068, the calibrator 
used was HD 10380 (or v Psc), a K3IIIb star whose angular 
diameter, derived from a photometric scale, is 2.99 ± 0.15 mas 
(Perrin et al. 1998). 

The interferometric response of the instrument, the so- 
called transfer function, was calculated from the calibrator ob¬ 
servations. The transfer function is the ratio of the observed 
fringe contrast by the expected visibility of the calibrator de¬ 
duced from the estimated diameter. Plots of the transfer func¬ 
tion as a function of wavelength are also informative about 

1 This software is now available through the JMMC website: 
http ://mariotti .uj f-grenoble. fr/ 
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Fig. 1. Visibilities for the four projected baselines as a function of wavelengths. Solid lines correspond to MIDI measurements 
while dashed lines are related to fits by the radiative transfer between two spherical components model (see sect.[3]an dg). 


the quality of the data with respect to turbulent conditions 
and background subtraction. Since sensitivity to turbulence de¬ 
creases with wavelength, the transfer function increases with 
wavelength. Excessive turbulence or misalignments cause a 
large gradient between 8 and 13 p m. Such data are rejected. 

The observations of the source are generally bracketed by 
observations of the calibrator star. The source fringe contrasts 
are obtained by normalizing the fringe contrast by the transfer 
function interpolated at the time the source was observed. A se¬ 
lection is then applied on the histograms of the squared fringe 
contrasts of NGC 1068 and of the calibration star. Sequences of 
fringe contrasts should have a rough gaussian distribution un¬ 
der steady turbulence conditions. Sequences of data for which 
histograms were departing too clearly from gaussian shape 
were rejected. 

Actually, two calibration stars were observed, T 4 Eri and 
HD 10380. A detailed look at the data led us not to use transfer 
functions derived from T 4 Eri since they strongly increase up 
to 11 pm and then saturate. As it was the only calibrator star 
used during the night of November 10th (except for a single se¬ 
ries), and as data obtained on Betelgeuse were also rejected for 
this same night (see Perrin et al. 12006k we did not use the data 
of November 10th for our present study of NGC 1068. In the 
same way, we were led to reject the data of June 15th since they 
most probably have been affected by large turbulence effects. A 
detailed look at the squared fringe contrast histograms also led 
us to reject some data from the two remaining nights. Table Q 
and Table|2]summarize the observations eventually kept for the 
present study. The combination of adaptive optics and spatial 
filtering will eventually solve this issues and improve the qual¬ 
ity of the data. 

Visibility measurements of NGC 1068 have been averaged 
by bins of projected baselines. A consistency check has been 


performed in each bin to assess the quality of error bars. Ad- 
hoc error bars have been computed to be consistent with the 
dispersion of visibilities when error bars were initially too 
small. The four bins are the visibility points used in the fol¬ 
lowing and are presented in Fig|7| The final absolute accuracy 
on visibility estimates is on the order of 3 percent. Our pro¬ 
cessing of the data leads to a more extended data set than that 
of Jaffe et al. ( 2004 ) who kept two visibility points at baselines 
of 42 and 72 m. 

Two main features are striking on Figin First, the visi¬ 
bilities are far smaller than 1, which means that NGC 1068 
is well resolved by MIDI. Second, visibilities are obviously 
wavelength dependent, a property which has to be taken into 
account in the modeling of the source. 

2 . 2 . Spectral Energy Distribution calibration 

The Spectral Energy Distribution (SED) as seen by MIDI is of 
prime importance to model the spectral and spatial properties 
of the source. The SED of NGC 1068 is derived from the chop¬ 
ping sequences of the source and of the calibrator. However, 
no observed SED is available for HD 10380 in the literature. 
Therefore, we made use of the calibrated one computed by 
Cohen et al. The integrated flux has been normalized to 

that of HD 10380. The two SEDs are consistent to within 7%. 
We eventually used the mean of these two spectra as a SED for 
HD 10380. The difference between the two SEDs was kept as 
the error on the final SED of HD 10380. A spectrophotometric 
response of MIDI was calibrated on HD 10380 using this SED. 

Two spectra of NGC 1068 were independently derived 
from photometries measured on June 16th and November 9th, 
2003. Error bars on photometry measurements are set from the 
difference between the spectra measured with the two UTs. The 









4 


Poncelet et al.: A new analysis of the nucleus of NGC 1068 with MIDI observations 


final MIDI SED of NGC 1068 is derived from the mean of the 
spectra of June and November and from the spectrophotometric 
response. The difference between the calibrated SEDs is used 
as an error estimate. 

As shown in Fig. |3 the photometric measurements of 
NGC 1068 and HD 10380 are strongly wavelength-dependent. 
The characteristic loss of signal around 9.5 fim is due to ab¬ 
sorption by atmospheric ozone. To avoid too large noise sen¬ 
sitivity, we do not take into account regions where photometry 
measurements are smaller than half the maximum value. These 
regions are greyed on our final spectrum of NGC 1068. The 
SEDs of June, November and the average one are presented on 

Fig E] 



Wavelength (microns) 

Fig. 2. NGC 1068 photometry measurements obtained with 
MIDI on June and November 2003. Greyed areas are related 
to wavelength ranges where photometries are smaller than half 
of their maximum value. 


The SED of NGC 1068 we derived from the MIDI data is 
different from the one presented by Jaffe et al. (ESS: quan¬ 
titatively, the former is in the range 13 to 30 Jy whereas the 
latter is rather between 8 and 15 Jy. Moreover the slope of the 
SED we derived after 11 yum is steeper than the one of Jaffe et 
al. Besides, we have validated our calibration procedure on an 
independent data set acquired on Betelgeuse. We have checked 
that we could reproduce the spectrum of Verhoelst et al. (2005) 
to an excellent accuracy. Besides, we have compared our spec¬ 
trum to the one available from the ISO database (see Fig. 0. 
The field of view of ISO is much larger than that of MIDI and 
ISO is detecting the dust shell whose inner edge starts at 0.5” 
(Danchi et al. 1994). The emission feature of silicate dust at 
9.7 /im is therefore present in the ISO spectrum but not in the 
MIDI spectrum whose field of view is reduced to less than 0.3”. 
Apart from dust, the continua should however agree. We find a 
maximum difference of 16% between our calibrated SED and 
the ISO spectrum outside the silicate dust feature which makes 
us confident in our calibration. 



Fig. 4. Comparison between the SEDs of Betelgeuse derived 
from MIDI photometry measurements (thick solid line) and 
the ISO-SWS spectrum (thin dotted line). This result on 
Betelgeuse makes us confident with our calibration procedure 
and with the fluxes we derive for NGC 1068 (see text in 
sect. 12 .2> . 

3. The different steps of data modeling 

As only a few MIDI data points on NGC 1068 are avail¬ 
able so far, models were computed from a rigorous stepwise 
approach on the data rather than from a priori knowledges 
about the source. Models remain also simple to avoid over¬ 
interpretations and to keep a reasonnable number of degrees 
of freedom. 

The characteristic apparent size of the object was first 
gauged by fitting visibilities with a uniform disk model whose 
diameter is wavelength dependent. The 39 m baseline data 
on the one hand and the 45-46 m baseline data on the other 
hand were fitted independently since these two sets of projected 
baselines are separated by 10° in azimuth. The 78 m data were 
not used at this stage since they would lay in the second lobe 
of the uniform disk visibility function, which model has a poor 
chance to be the right one at high spatial frequency. Therefore, 
these data would bias the measurement of the characteristic size 
of the source, which is the main goal of this first procedure. 

Fig.EJclearly shows the wavelength dependency of the ap¬ 
parent size of the object. One may notice that the curves re¬ 
lated to the two different baselines do not overlap with a dif¬ 
ference ranging from 2 to 7 mas. First, this may be due to the 
model which is clearly too simple. This discrepancy could then 
be interpreted as an asymmetry of the nucleus of NGC 1068, 
which one is well known at larger angular scales. Thus, to pos¬ 
sibly highlight some elongation and orientation of the source, 
an attempt was made to fit all the visibility points at the same 
time with a uniform elliptical disk. This led to a degeneracy 
of parameters which means that extracting informations on the 
shape and orientation of the source is not possible with the 
present MIDI data set. As a consequence, geometries used in 
the following modelings were restricted to spherical symmetry. 
The two typical sizes of ~ 30 mas and ~ 55-60 mas at 8 yum and 
after 9.5 yum respectively shown up in Eig.0suggest to elabo¬ 
rate a further model by the addition of a second uniform disk. 
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Fig. 3. SEDs of NGC 1068 derived from MIDI photometry measurements. Black solid line and vertical error barrs are related 
to the resulting SED of NGC 1068, which is the mean between those of June and November. The red solid line represent the 
fit of the SED obtained with the model of radiative transfer between two spherical components (see sect. EJandEJ- Greyed area 
correspond to wavelength ranges where photometry measurements are weak (see fig. El- 



Wavelength (microns) 



Fig. 5. Evolution of the diameter derived from the uniform disk 
model for two different configurations of the projected baseline 
(up: 39.59m alone; down: 45.54, 46.63 and 45.7m together). 


Fig. 6. Evolution of the diameters of the inner disk (solid line) 
and the extended disk (dashed line) derived from the two uni¬ 
form disks model. 


This model counts now three free parameters per wavelength, 
being the diameters of each disk and the flux ratio between 
them. Fits of visibilities are improved, which strengthens the 
validity of the two-disk model. Diameter variations with wave¬ 
length are plotted on Fig. [6] 

From these previous results, the third step in the modeling 
is to consider the radiative transfer between the two spherical 
components, namely a compact inner core (corresponding to 
the inner disk of diameter 0 ln ) which radiates as a black body 
with temperature 7’ ln and a layer (described by the outer disk 


of diameter 0i ayer ) around the compact core which radiates as a 
thermal grey-body at 7) avcr . The absorption in the layer and its 
emission are accounted for by its optical depth r(/l). No scat¬ 
tering is considered here. For the sake of simplicity, we used 
a model in which the layer is geometrically thin. This model 
was successfully applied to Mira stars and red supergiants in 
the near-infrared (see e.g. Perrin et al. 120041 . The interest of 
this model is that the radiative transfer can be calculated ana¬ 
lytically and that there is no major difference with a model for 
which the base of the layer would be at the edge of the compact 
inner core. 
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Thus, if 0 j n 3; 0i a yer sin# (where 6 is the angle between the 
line of sight and the radius vector of the two spheres), both 
contributions of the compact core and of the layer are observed. 
Then the emerging intensity can be expressed as: 

I(A,d) ~ B(A, T in ) e~ ™ )/coifl 
+ B(A , T layer ) [1 - e- r( ' l)/a " fl ] 

In the other case. 


However, several local minima appear in the x 2 - 
hypersurface in temperature space, leading to a degeneracy of 
these parameters. Indeed, visibilities are not sensitive to the ab¬ 
solute flux of each component, but to the relative flux between 
them. 

The degeneracy on temperatures is overcome using the 
MIDI SED of NGC 1068. The criterion estimating the good¬ 
ness of a fit of the SED, with a given couple of temperatures, is 
given by: 


I(A,0) ~ B(A, T layer ) [1 - e~ 2 l -W)/«™ s ] 

where B(/l,r) is the Planck function. 


X Spec 


I 


SED obs (i) -SED bt (i) 


0"sed(O 


Thus, as the model was chosen to remain rather simple, 
the number of free parameters was limited to five at each 
wavelength. Diameters and temperatures are global parameters, 
while the optical depth, chosen to be wavelength dependent, 
allows to account for the visibility variations and for the SED. 
The total number of visibility data is 4xV„i while the number of 
free parameters is then 4 + N,\ (N,i = 109 being the number of 
spectral channels). The diameters are well constrained by the 
visibilities while the latter only provide a one-to-one relation 
between the temperatures. Compatibility with the SED allows 
to determine a unique couple of temperatures. 

The criterion to be minimized is the^ 2 , written as: 

N 

T„, = £ 

;=i 

where M(0 in , 0] ayer , T m , T\ ayeT , tx, St) is the modeled 
squared visibility function, which also depends on the spatial 
frequency .S', given by: .S', = B/A, (B being the projected 
baseline). 

There are several local minima on the x 2 hyper-surface in 
parameter space. Therefore the search for optimal parameter 
values has to be carried out in several steps. To reach the solu¬ 
tion diameters, the first step consists in the exploration of a 2-D 
grid in ( 0 j n ; 0 i ayer ) space. For each point on this grid, the mini¬ 
mization is performed on the temperatures and on a global opti¬ 
cal depth (in other words, a mean optical depth non-dependent 
on wavelength at this point of the minimization). For the sec¬ 
ond step, the same procedure is applied but on a 2-D grid in 
(Tinflayer) space. The third and last step consists in running the 
minimization on optical depth with respect to wavelength, set¬ 
ting diameters and temperatures at the optimum values found 
previously. Ideally, steps should have been iterated but this did 
not turn out to be necessary. 


Vi M( 0 in, 0layer> ^layer? S j) 

0~i 


4. Results on the NGC 1068 mid-IR core 

Visibility measurements actually well constrain the geomet¬ 
rical parameters since the first step of the x 2 minimization 
rapidly converges to the optimum values of geometrical param¬ 
eters, which are: 

0 ; n = 35.3 ± 3.7 mas 
01ayer = 83.1 +6.1 mas 


The optimal couple of temperatures then obtained is: 

Tin = 361 ± 12 K 

Tjayer = 226 + 8 K 

(bvis ^Reduced = ^3.9 

V 2 s P ec - 16.4 

Xy ls and Aspec can be minimized independently as long as 
minimization of the second one is performed over the local 
minima of the first one. 

To estimate the confidence interval on parameters, we treat 
the two couples of global parameters (i.e. diameters and tem¬ 
peratures) independently (e.g. Bevington & Robinson 
For geometrical parameters, we look at the joint confidence in¬ 
terval given by a variation of 1 with respect to^y is /( / \'y is ) m i n , as 
a function of 0 i n and 0 i ayer , keeping the other parameters fixed 
at their optimum values. This is equivalent to increase the error 

bars on visibilities by Jx 2 / Ofvis )™ 11 an d to reaccount for the 
distance to the model. 

For temperatures, the confidence interval corresponds to the 
intersection between the ^ 2 |s trough in f/) n ;7'| ayc . r ) space and the 
joint confidence interval on temperatures which give a reason¬ 
able fit of the SED. In order to obtain this last joint interval, 
we look at a variation of 1 on the value of Tspec/Ofspec)™ n as 
a function of T m and ri ayer . Values of x 2 relating to the fits of 
the visibilities and of the spectrum are quite high, leading to 
an under-estimation of the confidence interval on parameters. 
Nevertheless, this must be due on the one hand to the simplic¬ 
ity of the model considered and on the other hand to the fits of 
the 45 and 46 m visbility points between 8 and 8.5 fjm which 
are poor compared to eslewhere in the A'-band (see Fig. |TJ. 
Also note on Fig.^that the modeled visibility curves have the 
same kind of wavelength dependent fluctuations as the mea¬ 
surements. This is due to the third step of the x 1 minimisation 
which improves the fit by converging to the optimal value of 
the optical depth at each wavelength. 

Fig m presents several fits of visibility measurements as 
a function of the projected baseline at different fixed wave¬ 
lengths. The growing influence of the extended layer with in¬ 
creasing wavelength is highlighted on the one hand by the flat¬ 
ness of the visibility function at high spatial frequencies and 
on the other hand by the strong slope at low spatial frequen¬ 
cies. Moreover, the figures underline the embarrassing lack of 
visibility points at low projected baselines. As the slope at low 
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Fig. 7. Fit of the visibilities as a function of the projected baseline and at given wavelengths with the radiative transfer between 
two spherical component model. 


spatial frequencies is really model-dependent, these points are 
needed to put stronger constrains on the modeling. 

The fit of the MIDI spectrum is presented on Fig. 0 This 
shows that not only the simple radiative transfer model ac¬ 
counts for MIDI visibility measurements, but also it gives a 
good representation of the SED across the full A'-hand. Fig. [8] 
presents the contribution of each component to the spectrum. 
We observe that the inner compact component is the only one 
to contribute to the SED until 8.5 fim while the flux emitted 
by the layer begins to dominate above 9.4 fin l. This figure also 
shows that, in the framework of our model, the gap of the SED 
around 10 fim is not due to the crossing of the two black body 
emissions, but indeed to a peak in the absorption by the layer. 

Therefore, the originality of the present study resides in its 
capability to provide the evolution of the optical depth of the 
extended component as a function of wavelength in the /V-band, 
without any a priori assumption on it like its variation law or 
the composition of the extended layer. The optical depth result¬ 
ing from the whole analysis is plotted as a function of wave¬ 
length in Fig.|b](dark solid line) and will be discussed further 
in sect. 15.31 


5. Discussion 

5.1. Spatial description of the source 

As the number of visilitity measurements is extremely small 
and as the (u,v) plane is not sufficiently covered, models able to 
fit the MIDI data points are degenerated. This lack of sensitivity 
in the models highlights the importance of following a stepwise 
approach from the data. 

The MIDI visibility points can be fitted by a different and 
a priori model that consists in a two gaussian elliptical disk 
model (Jaffe et al. 120041 . The derived sizes are ~ 10 x 12 mas 
for an inner hot, optically thick component and ~ 30 x 49 mas 
for a warm component. The factor of 2 between these values 
and the ones we derived with our model is actually explained by 
the differences between the respective geometries considered. 

Comparing mid- and near-IR observations, it appears that 
the largest A'-hand structure corresponds in size to the small¬ 
est structure detected in the mid-IR with MIDI. From A'-hand 
bispectrum speckle interferometry, Wittkowski et al. ( 19981 in¬ 
terpreted the observed central core of NGC 1068 as a slightly 
resolved gaussian disk of FWHM ~ 30 mas. Speckle imaging 
from the Keck Observatory also showed that the flux emitted 
by the core of the AGN may be accounted for by an unre¬ 
solved point source whose size is smaller than ~ 20 mas and 
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Fig. 8. Contributions of each spherical components of the ra¬ 
diative transfer model to the SED. Red thin dotted line-, flux 
emitted by the non-absorbed inner component. Blue thick dot¬ 
ted line: fraction of the flux emitted by the inner component and 
absorbed through the layer. Red thin dashed line : flux emitted 
by the layer without its self-absorption considered. Green thick 
dashed line : flux emitted by the layer and self-absorbed by it. 


which is surrounded by an emitting region extending up to 
10 pc (Weinberger et al. Wittkowski et al. (1200411 com¬ 

bined near-IR speckle data with the K -band visibility point ob¬ 
tained with VINCI and interpreted the whole dataset with a 
two-component model: a small one having a size of ~ 5 mas (or 
0.4 pc) responsible for 40% of the X-band flux and an extended 
one of ~ 40 mas. In addition, near-IR bispectrum speckle inter¬ 
ferometric observations led to a N-W extended structure in K'- 
band having a size of 39 x 18 mas and a RA. ~ -16° (Weigelt et 
al. 2004k With regard to the AO observations performed with 
NACO at VLT, the k'-band core is slightly N-S elongated with a 
FWHM of 67 mas (Rouan et al.. 12004 ). The spatial description 
of the source derived from speckle observations is also fully 
consistent with the best resolution AO observations performed 
with NAOS-CONICA which show a resolved structure in the 
/Cv-band, elongated along the same P.A. ~ -16° and of FWHM 
~ 30x <15 mas after deconvolution (80 x 65 mas before de- 
convolution; Gratadour et al. 12005k The sizes derived from K- 
band and ,V-band observations underline the consistency be¬ 
tween the near-IR and the mid-IR observations of the nucleus 
of NGC 1068. 

According to Maloney the dusty torus has suit¬ 

able temperatures (T ~ 250 K) and FF molecules densities 
(«h 2 ~ 10 8 - 10 10 cm -3 ), in order to produce the 6ie —» 523 
rotational transition at 22 GHz of the H 2 O masers (e.g., Elitzur 
119921 chap. 10). H 2 O masers can be produced by the absorp¬ 
tion by a H 2 O molecule of an X-ray emitted by the central en¬ 
gine of the AGN, and the 22 GHz radio emission is expected 
to trace the inner edge of the torus. Sub-milliarcsecond angu¬ 
lar resolution observations made on NGC 1068 with the Very 
Fong Baseline Interferometer (VFBI) showed that the red- 
shifted emission is divided in four groups, linearly distributed 
between 6 and 14 mas from the core along a direction at 45° 
from the radio synchrotron jet axis (Greenhill et al. fl996t . 


Additional observations of the blueshifted emission showed 
that masers are actually spatially distributed on a thin ring in¬ 
clined at -45°, on scales extending from 0.65 to 1 pc (i.e., 9.3 
and 15.7 mas; Greenhill & Gwinn fHl97l i. Thus, the ring traced 
by the H 2 O masers seems to correspond to the inner edge of 
the dusty torus (Kartje et al. 11999k This is actually what we 
are observing since the masers are located at the same place 
as the V-band structure derived from the VINCI/Speckle ob¬ 
servations and also observed with NAOS-CONICA, and just 
inside the extended component derived from our modeling of 
the MIDI data. 

Moreover, as claimed by Rouan et al. 1 2004 ) and Gratadour 
et al. (.2 005 1 . and according to current estimates of the dust sub¬ 
limation radius, graphite grains are able to survive at such small 
distances from the core of the nucleus of NGC 1068, so that this 
dust species could be part of the observed V-band structure. 
Indeed, a theoretical estimate of the dust sublimation radius is 
given by: 

r sub — 1-3 (Cuv,46) 1/2 (Tisoo) 2 8 PC (1) 

where Cuv ,46 is the UV luminosity (in units of 10 46 erg.s -1 ) 
and Tigoo is the sublimation temperature of the dust species 
considered (in units of 1500 K; Barvainis 119871 . Upper lim¬ 
its are then obtained by setting the UV luminosity of the cen¬ 
tral source to the observed bolometric luminosity of the nu¬ 
cleus and by considering sublimation temperatures of 1400 K 
and 1750 K for silicates and graphite grains respectively. Pier 
et al. ( 1994 ) investigated the bolometric luminosity of the nu¬ 
cleus of NGC 1068. Depending on the fraction of nuclear flux 
reflected into our line of sight which has been determined 
by several different ways, the bolometric luminosity of the 
nucleus of NGC 1068 ranges from ~ 6.8 x 10 43 erg.s -i to 
~ 3.4 x 10 45 erg.s -1 , the most probable value being on the order 
of ~ 3.4x 10 44 erg.s -1 (Pier et al. 119941 ). These dispersed values 
lead to a sublimation radius ranging from ~ 0.07 to 0.50 pc (i.e. 
~ 1 to 7 mas) for graphite grains and from ~ 0.13 to 0.91 pc (i.e. 
~ 1.8 to 12.6 mas) for silicate grains, the most probable values 
being ~ 0.16 pc and ~ 0.3 pc (or ~ 2.2 mas and ~ 4.2 mas) 
respectively. 

As a conclusion, in spite of the great uncertainties on the es¬ 
timation of sublimation radii, it is remarkable that H 2 O masers 
seem to be confined in an area located between the sublimation 
radius of silicate grains and the outer edge of the inner compo¬ 
nent we have inferred from mid-IR data. In addition, according 
to near-IR observations of the nucleus of NGC 1068, it seems 
that the dust surrounding the central engine is divided in several 
layers of which the first one - emitting in X-band and probably 
extending from the sublimation radius of graphite grains up to 
about ~ 1.2 pc - could be mainly composed of graphite grains 
and would harbour the observed H 2 O masers. 

5.2. Temperature of the MIDI compact component 

The temperatures we derive are much cooler than those of Jaffe 
et al. mm . There are different reasons to these discrepan¬ 
cies. The most important one comes from the difference in 
the geometries considered that leads to a factor of 2 between 
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the respective sizes derived and consequently to a factor of al¬ 
most 0.5 on the temperatures. Moreover, our spectral calibra¬ 
tion yields twice as much flux. Finally, contrary to Jaffe et al., 
we do not consider absorption in the inner compact compo¬ 
nent which we modeled as a pure black body. Indeed, for Jaffe 
et al., the optical depth associated to silicates in the inner hot 
component is about 2.1, meaning that the radiation emitted by 
the central engine is greatly attenuated in this optically thick 
medium. In this paper, we do not consider this additive param¬ 
eter and the opacity of the inner component affects the final 
temperature. 

The temperatures we derive from the radiative transfer 
modeling seem very low with respect to dust sublimation tem¬ 
peratures commonly assumed for graphite grains and silicates 
which are on the order of 1700 and 1450 K respectively. 
However, following the approach of Krolik ( 1999 ;) who, for 
first rough estimates, assumes thermal emissions by spherical 
black bodies, the typical angular size of each component inside 
an AGN is given by the expression: 


0 = 87 


Fobs 


io- 


x erg.cm 2 .s 


1/2 


F e ff 

Wk 


x-2 


p as 


( 2 ) 


Hence, for typical sizes of 35 and 83 mas and a flux of 
~ 17 Jy at 8 pm and ~ 30 Jy at 13 pm (according to the 
MIDI spectrum of NGC 1068), the typical temperatures ob¬ 
tained with this formula are ~ 260 and 160 K. These values 
are consistent with those we obtain with the radiative trans¬ 
fer model. This shows that considering this simple description, 
even at small distance from the central engine, temperature falls 
quickly. Thus, the apparent weakness of temperatures derived 
from our modeling is finally not surprising. 

Furthermore, according to SED measurements of the cen¬ 
tral elliptical region of 290 x 180 mas of the nucleus of 
NGC 1068 carried out at the Subaru telescopes with a 0.1” 
spatial resolution, Tomono et al. ( 200 11) derived a tempera¬ 
ture of ~ 234 K and an optical depth ~ 0.92, from a 

model of modified greybody radiation absorbed by silicates. 
Temperature for such a large area is consistent with those we 
obtain and the optical depth we are able to derive at 9.7 pm is 
fairly close. 


5.3. Dust absorption in the layer 

The optical depth as a function of wavelength across the N- 
band is deduced from the optimisation of the fit of visibility 
measurements. Therefore, it is rather model-dependent. 

However, a detailed study of the shape of the optical depth 
obtained with our model is interesting. Indeed, looking at 
Fig- 0 it seems that the strong slope between 8.5 and 9.5 pm 
and the bump around 10 pm can be the signature of the pres¬ 
ence of amorphous silicates in the composition of the layer 
(in comparison with the absorption features at 10 pm associ¬ 
ated with the Sgr A* region observed with ISO for instance; 
Kemper et al. * 2004 ). Moreover, given the sublimation radius 
of silicates which is on the order of 0.64 pc (or ~ 8.9 mas; 
see sect. 15. 1> . this type of dust is actually able to survive in the 
layer. Therefore, to verify the probable presence of silicates, we 


tried to fit the shape of the derived optical depth with DUSTY, 
a radiative transfer code into a dusty layer, which includes ab¬ 
sorption, emission and scattering by dust (Ivezic, Nenkova & 
Elitzur, This code handles spherical geometry and has 

been used by Imanishi & Ueno (2 000 ) to model the Cygnus A 
nucleus absorption feature at 9.7 pm. It is important to under¬ 
line here that if the torus is viewed almost edge-on (the com¬ 
mon assumption for NGC 1068), the use of spherical symmetry 
in the modeling does not affect the derived optical depth since 
the thickness of material viewed along the line of sight remains 
the same as in the case of a torus-like geometry. 

Next, Maiolino et al. ( 2001a 2001b ) claim that dust in the 
circumnuclear region of AGNs must have different properties 
than in the galactic interstellar medium (ISM) and that sili¬ 
cate grain sizes must be larger than 3 pm in order to explain 
the lack of silicate features in ISO spectra of Seyfert II nu¬ 
clei. This is confirmed by the computation of the T3.4/T9.7 ratio 
which is ~ 0.11 for NGC 1068 (Dartois et al. 2 0041 ). Indeed, 
this ratio would be ~ 0.06-0.07 assuming the Av/rgj ratio of 
both the galactic diffuse interstellar medium (DISM) in the lo¬ 
cal 3 kpc around the Sun (Roche & Aitkcn 119841 and in the 
direction of the Galactic Centre (Roche & Aitken 1X9851 . and 
the Av/t3.4 ratio in these different directions (Pendleton et al. 
Il994t . In addition to grain sizes, such a difference in extinction 
laws of the Galactic ISM and Seyfert II AGNs can also result 
from a temperature gradient in the obscuring dust (Imanishi 
12000) . Evidences of the flattening of the silicates optical depth 
with the growth of grain sizes are presented by Bouwman et al. 
( 2001 ) and Kemper et al. ( [2004) . 

Following these considerations, tentative fits of the de¬ 
rived optical depth with several simple compositions of amor¬ 
phous silicates and various grain sizes provided with the code 
DUSTY are presented in Fig.loland l 1 01 We are considering the 
MRN distribution (Mathis, Rumpl & Nordsieck ll977l i which 
consists in a power law (index = -3.5) distribution of grain 
sizes ranging from a minimum value of 0.05 micron to vari¬ 
ous maximum sizes. Species of dust available in DUSTY are 
“astronomical silicates” and graphite grains from Draine & 
Lee 119841 , referred as Sil-DL and Grf-DL thereafter), warm 
oxygen-deficient and cool oxygen-rich silicates - representa¬ 
tive of circumstellar and interstellar silicates - from Ossenkopf 
( 1992 1 referred as Sil-Ow and Sil-Oc), a Silicon carbid from 
Pegourie ( 19881 referred as SiC-Pg), and amorphous carbon 
(referred as amC). A first remark is that compositions of Sil-Ow 
or Sil-Oc alone will provide the same optical depth (in shape 
and amount); the addition of grf-DL or amC in the composition 
of the dust has no effect on the shape of the optical depth law 
since absorption or emission features of these species appear 
around 4 pm. We therefore only considered Sil-Ow, Sil-DL and 
SiC-Pg compounds thereafter. 

Sil-DL alone provides a good fit of the optical depth after 
10 pm while the 10 pm silicate feature is globally shifted to¬ 
wards lower wavelengths (see Fig |9}. On the contrary, when 
SiC-Pg are added, the bump is shifted to the right and it does 
not account for the shape around 11 pm. Best fits are obtained 
with a mixture of Sil-Ow with a small amount of Sil-DL. 

Next, for the study of grain sizes, we considered two of 
the well fitting compositions: 50% Sil-DL + 50% Sil-Ow and 
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Fig. 9. Optical depth as a function of the wavelength derived 
from the radiative transfer model (thick line); fits provided by 
the DUSTY code with different compositions of silicates are 
superposed. Sil-Ow refers to warm siliactes from Ossenkopf 
G23; Sil-DL refers to silicates from Drain & Lee Q2E3; 
SiC-Pg refers to a silicon carbide from Pegourie ([7988). 

100% Sil-Ow. We first considered a single grain size distri¬ 
bution. As shown by Kemper at al. CTS3I . we observed that 
the silicate feature is significantly shifted towards longer wave¬ 
lengths as the grain size increases. Best fits are provided by 
a size of 1 micron. Nevertheless, the MRN distribution pro¬ 
vides better fits. According to Fig. rroiwhich presents compar¬ 
isons of fits provided by various maximum grain sizes (rang¬ 
ing from 0.25 to 10 microns), distributions of relatively small 
grains (with sizes smaller than 3 microns) account better for 
the bump around 10 pm while larger grains flatten the optical 
depth law beyond 10 pm. 

The quality of the fit of the optical depth is poorer at the 
edges of the band, and the optical depth derived from the mod¬ 
eling appears to be more asymmetric than what would be ex¬ 
pected from DUSTY. There are several reasons for that: first, 
the model considered here for NGC 1068 is a fairly simple 
one. Consequently it is not meant to be a real description of the 
source compared to its a priori complexity. Thus, differences 
between the real geometry and physical properties of the source 
and the model will result in inaccuracies in the estimate of the 
optical depth. Second, differences may arise from the physical 
properties which are modeled. No scattering is taken into ac¬ 
count in the two component radiative transfer model used to fit 
the MIDI data. Moreover, following Kemper et al. ( 1200 4) who 
assume that the spectra of optical depth are all the narrower as 
emission by silicates increases, the narrowness of the shape of 
the optical depth provided by DUSTY beyond 11.5 pm could 
highlight that our model under-estimate silicate emissions in 
the layer. 

However, if we now compare the evolution of optical depth 
with the one provided by the ISO observations of Sgr A* 
(Kemper et al. 120041 see Fig. m, it appears that, although 
the silicates peak is displaced with respect to wavelength, both 
plots have similar slopes at low wavelengths. Besides, despite 
there are no evidences yet of the presence of PAH (Polycyclic 



Fig. 10. Fit of the optical depth provided by the DUSTY code 
for various maximum sizes of grains, for a MRN (Mathis, 
Rumpl & Nordsieck ll977t grain size distribution and two dif¬ 
ferent compositions. Top: 50% silicates from Ossenkopf ( 1992 ) 
and 50% silicates from Drain & Lee am and superposi¬ 
tion with silicate features intrinsic to the Sgr A* region. Down: 
100% silicates from Ossenkopf dl992lb . 

Aromatic Hydrocarbon molecules) features in the ISO spec¬ 
trum of NGC 1068 (at 7.7, 8.6 or 11.3 pm; Sturm et al.EEP, 
this trail should be investigated in a further study. Similarly, the 
strong pure-rotation lines of H 2 O could significantly contribute 
to opacity across the /V-band and could also provide an expla¬ 
nation for the wing in the upper half of the band. Therefore, 
as these molecules would be able to survive in an environ¬ 
ment such as the dusty layer, their presence could also affect 
the shape of the optical depth and be responsible for the non 
perfect fit. 

Nevertheless, at this point conclusions can already be 
drawn from the present study of the optical depth. The main 
features of the dust around the core of NGC 1068 could be: 

- a mixture of amorphous silicates from Draine and Lee 
dT984il with a large fraction of cosmic silicates from 
Ossenkopf 11 992 1). 

- small grains, with sizes ranging from 0.05 to 3 microns. 

The properties derived from the optical depth deduced from the 
MIDI observations show that dust around the central engine of 
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Fig. 11. Comparison between the results derived from the ra¬ 
diative transfer model between two spherical components ap¬ 
plied to MIDI data with others observations. 


the nucleus of NGC 1068 has been detected. This is the first 
direct evidence of the presence of amorphous silicates in the 
composition of the dusty environment at radial distances of 1 
to 2 pc from the central engine. 

5.4. Proposed description of the dust distribution 

From the above discussion, a picture of the distribution of dust 
inside the core of NGC 1068 can be drawn (see the scheme 
in Fig. IIH . It is a layered structured. One of the layers is de¬ 
tected in the K- and /V-bands and extends from the sublimation 
radius of graphite grains up to about 1.2 pc. It could be made 
of graphite grains, and silicates in its outer parts. The second 
one emitting in the /V-band and extending from 1.2 pc up to 
about 2.7 pc, is well resolved by MIDI and seems to be warmer 
(about 200 K) and made of amorphous silicates. This picture is 
consistent with IDO maser emission since masers need a dusty 
environment in order to survive in the hard U V field radiated by 
the central engine. The first layer would then act as a protecting 
shield and would be a favorable environment for water maser 
emissions. It also provides an explanation for the low temper¬ 
atures derived from the radiative transfer model applied to the 
MIDI observations. Indeed, there would be a strong radial tem¬ 
perature gradient inside the inner layer. 


6. Conclusions 

We present an independent analysis of the first VLTI MIDI 
observations of the nucleus of NGC 1068 obtained in 2003. 
The small visibilities measured across the /V-band show that 
the core of NGC 1068 is well resolved in the mid-IR. 

After a severe selection of the data, different basic models 
with an increasing number of free parameters have been tried 
to fit both interferometric and spectral data. As a final scenario, 
we adopt a simple radiative transfer model with two concen¬ 
tric spherical components. This model accounts well for the 
MIDI visibilities and SED across the /V-band. The inner com¬ 
pact component has a radius of about 17 + 2 mas and a tem¬ 
perature of about 361 ± 12 K, and the outer layer extends up 
to 41 ±3 mas from the core of the nucleus and has an average 
temperature of about 226 + 8 K. Results are consistent with 
observations in other spectral bands of the thermal IR and with 
FRO masers emissions. One specificity of this approach is that 
it provides the evolution of the layer optical depth across the N 
band, without requiring assumptions about its nature. Although 
this parameter is rather model-dependent, its shape seems to be 
characteristic of amorphous silicate grains. Thus, the conclu¬ 
sion to be drawn from the analysis of the data with this simple 
model is that MIDI has actually observed the dust distribution 
around the core of the nucleus of NGC 1068 emitting in the 
A'-band. 

Besides, this study highlights the degeneracy in the num¬ 
bers of models being able to describe the current MIDI ob¬ 
servations. This should be overcome thanks to various coming 
near- and mid-IR observations at the VLT. Concerning MIDI, 
the coverage of the (u,v) plane needs to be increased in order 
to get visibility points for several new lengths and orientations 
of the projected baseline. Observations with MIDI should be¬ 
sides put additional constrains on the modeling and lead to a 
convergent description of the morphology and of the physical 
properties of the dust distribution. Next, as the slope of the vis¬ 
ibility function considered in the model presented by Jaffe et 
al. 12004 > differs from ours at low spatial frequency, visibility 
points at low projected baseline should lead to the choice of the 
most favorable model. As the shortest projected baseline reach¬ 
able with the VLTI and the UTs is ~ 30 m, one would have to 
make use of VISIR (the VLT imager and spectrometer for the 
mid-IR) to get 8.2 m baseline visibility points. 

This study is a very first step in the understanding of the 
dusty environment of the core of an AGN. It will be mandatory 
to test the approach presented here with observations of other 
AGNs such as Centaurus A or Circinus. 
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Table 1. Log of June 16, 2003 observation, performed between UT1 and UT3 


UT time 

Projected baseline 
(m) 

Baseline Azimuth 

(°) 

UT time 
(Calibrator 1) 

UT time 
(Calibrator 2) 

09:56:54.000 

78.6670 

0.036 

09:34:22.000 

- 

09:58:46.811 

78.6670 

0.036 

09:35:38.451 

- 

10:02:32.432 

78.6670 

0.036 

09:36:54.901 

- 

10:10:12.000 

78.6670 

0.036 



10:11:44.611 

78.7050 

0.036 



10:13:17.223 

78.7050 

0.036 



10:14:49.834 

78.7050 

0.036 




Table 2. Log of November 9, 2003 observations, performed between UT2 and UT3 


UT time 

Projected baseline 
(m) 

Baseline Azimuth 

(°) 

UT time 
(Calibrator 1) 

UT time 
(Calibrator 2) 


02:50:22.451 

39.5900 

35.16 

- 

03:59:13.432 


02:52:55.352 

39.5900 

35.16 

- 

04:00:29.882 





- 

04:01:46.333 


04:48:28.000 

45.5410 

44.60 

03:59:13.432 

05:17:54.000 


04:49:44.451 

45.5410 

44.60 


04:00:29.882 

05:19:10.450 

04:51:00.901 

45.5410 

44.60 


04:01:46.333 

05:20:26.901 

05:58:27.451 

46.6340 

45.90 

05:17:54.000 

06:33:28.000 


05:59:43.882 

46.6340 

45.90 

05:19:10.450 

06:34:44.450 





05:20:26.901 

06:36:00.901 


07:01:56.000 

45.7090 

44.76 

06:33:28.000 

_ 


07:03:12.450 

45.7090 

44.76 

06:34:44.450 

- 





06:36:00.901 

- 




